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Abstract 

Notch signaling is involved in cell fate choices during the embryonic development of Metazoa. Commonly, Notch signaling 
arises from the binding of the Notch receptor to its ligands in adjacent cells driving cell-to-cell communication. Yet, cell- 
autonomous control of Notch signaling through both ligand-dependent and ligand-independent mechanisms is known to 
occur as well. Examples include Notch signaling arising in the absence of ligand binding, and cis-inhibition of Notch 
signaling by titration of the Notch receptor upon binding to its ligands within a single cell. Increasing experimental 
evidences support that the binding of the Notch receptor with its ligands within a cell (cis-interactions) can also trigger a 
cell-autonomous Notch signal (cis-signaling), whose potential effects on cell fate decisions and patterning remain poorly 
understood. To address this question, herein we mathematically and computationally investigate the cell states arising from 
the combination of cis-signaling with additional Notch signaling sources, which are either cell-autonomous or involve cell- 
to-cell communication. Our study shows that cis-signaling can switch from driving cis-activation to effectively perform cis- 
inhibition and identifies under which conditions this switch occurs. This switch relies on the competition between Notch 
signaling sources, which share the same receptor but differ in their signaling efficiency. We propose that the role of cis- 
interactions and their signaling on fine-grained patterning and cell fate decisions is dependent on whether they drive cis- 
inhibition or cis-activation, which could be controlled during development. Specifically, cis-inhibition and not cis-activation 
facilitates patterning and enriches it by modulating the ratio of cells in the high-ligand expression state, by enabling 
additional periodic patterns like stripes and by allowing localized patterning highly sensitive to the precursor state and cell- 
autonomous bistability. Our study exemplifies the complexity of regulations when multiple signaling sources share the 
same receptor and provides the tools for their characterization. 
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Introduction 

The Notch signaling pathway mediates cell-to-cell communica- 
tion in several developmental contexts [1-4]. This communication 
occurs through the binding of the Notch receptor in a cell 
membrane to its ligand (e.g. Delta) in a neighboring cell [5], what 
herein we refer as trans-interactions. The bound complex is then 
cleaved and its intracellular domain (NICD) targets gene 
expression within the cell that harbors the bound receptor [6]. 
In the case of neural development, Notch signaling mediates 
lateral inhibition, which drives the selection of cells that express 
high levels of proneural genes and ultimately become neurons [7— 
11]. Notch signaling activity inhibits the proneural genes, which in 
turn activate Delta expression. As a result, a cell committed to the 
neural fate with high proneural gene and Delta levels inhibits its 
neighboring cells from adopting the same fate (i.e. it performs 
lateral inhibition) [12]. When all precursor cells are initially 
equivalent and signal similarly, mutual lateral inhibition arises and 
can drive a spontaneous spatially periodic selection of precursors 
by amplifying the initial small differences between them [13]. 



Different experimental evidences show that the Notch receptor 
can bind as well to the ligands when they are both in the same cell, 
what we refer as cis-interactions [5,14—18]. Cis-interactions drive a 
reduction of Notch signaling by sequestering the receptor and 
impeding its signaling [14,16,18-28]. This is known as cis- 
inhibition and its effects have started to be theoretically and 
computationally addressed too [25,29-38]. These studies have 
revealed that cis-inhibition can facilitate patterning by promoting 
faster responses, enhancing robustness and precision, and relaxing 
the constraints required for patterning [25,29,32,33,35,37,38]. 

Although cis-interactions commonly inhibit Notch signaling 
[14,16,18-28], Notch activity coming from cis-interactions has 
been proposed for specific scenarios [39—4-3]. For instance, 
Coumailleau et al. (2009) pointed to cell-autonomous Delta- 
dependent active Notch in Sara endosomes of Drosaphila bristle 
precursor cells. Guy et al. (2013) indicated that cell-autonomous 
Notch-mediated activation of the cell-cycle regulator c-Myc in 
mouse T-cells is impaired when ligand-receptor cis-binding is 
prevented. The mechanism by which cis-driven signaling can 
occur in these scenarios is still missing. According to the 
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mechanism proposed by Furfhauer and Gonzalez-Gaitan (2009), 
ligand-receptor binding within multivesicular endosomes could 
drive the release of Notch intracellular domain, driving Notch 
activation. This binding would occur in anti-parallel configura- 
tions (like for trans-interactions), as opposed to the parallel binding 
that is commonly associated to cis-interactions when occurring in 
the cell membrane and which is believed to prevent signaling. The 
interplay between Notch and the endocytic routes is starting to be 
uncovered and may shed light on this issue. 

In the present work we take advantage of computational and 
mathematical modeling to address the question of which would be 
the effect expected from cis-signaling when another signaling 
source that also uses the Notch receptor is acting (Fig. 1). Since 
both signaling sources use the Notch receptor, both can compete 
for it. The additional source of signaling, hereinafter referred to as 
primary signaling source, can be associated with trans-interactions. 
Yet, our modeling approach is not exclusive for such trans- 
interactions. The primary signal can be also driven by additional 
alternative mechanisms. In this context, recent work has shown 
that ligand-independent Notch activity can arise from the binding 
of Notch to other factors and from impaired endocytic regulation 
[39,40,44—54]. For instance, a recent study in Drosophila blood cells 
has detected a ligand-independent Notch signal that has a 
significant role in their development [52,53]. Our results show 
that when acting together with a primary signal, cis-signaling can 
act as cis-activation or as cis-inhibition. Competition between 
signaling sources underlies this switch. We establish under which 
conditions each regime arises. An extensive analysis of the 
parameter space shows that cis-inhibition enriches patterning, as 
opposed to cis-activation. Cis-inhibition promotes pattern multi- 
stability and modulates the selection of precursors. In addition, cis- 
inhibition facilitates cell-autonomous bistability. 

Results 

A simple model for lateral inhibition with cis-signaling 

We set a mathematical phenomenological model that includes 
two sources of Notch signaling: a primary signaling source and 
signaling driven by cis-interactions between the Notch receptor 
and its ligand within the same cell. The primary signal can be 
driven either by trans-interactions between the Notch receptor and 
its ligand in an adjacent cell (Fig. 1A) or through ligand- 
independent mechanisms (Fig. IB). We assumed that both the 
primary signaling source and cis-interactions drive the same type 
of signal. W e set primary Notch signaling to occur in a graded non 
switch-like fashion and to saturate to a maximal value, as recently 
experimentally reported for trans-interactions [25]. We assumed 
that cis-interactions can drive a graded increasing production of 
Notch signal activity with ligand up to saturation like the primary 
signaling source does (Methods). When both signaling sources are 
acting, the production of signal that each of them drives depends 
on the other source since both sources use the Notch receptor to 
signal. Accordingly, these productions in cell ( read (Methods): 

Pl '- = TT^Ta (1) 



where the primary source produces Notch signal activity at rate 




Figure 1. A model for Notch signaling driven by a primary 
signaling source and by cis-interactions. Cartoons of the Notch 
signaling components under study for (A) two adjacent cells that 
interact and for (B) an isolated cell. Black arrows stand for activation 
while red blunt arrows denote inhibition. (A) The ligand (red) in a cell 
binds the Notch receptor (blue) in a neighboring adjacent cell (trans- 
interactions). This elicits a Notch signal (NICD) that inhibits ligand 
production in the adjacent cell. The ligand can also bind the receptor 
within the same cell (cis-interactions) and drive Notch signaling at a 
different strength (dashed arrow). (B) A primary signaling source that is 
ligand-independent is depicted as well as signaling driven by cis- 
interactions. In both panels, the Notch signal inhibits the ligand 
through the proneural genes. 
doi:10.1371/journal.pone.0095744.g001 

P\i whereas cis-interactions produce it at rate Pii- In the above 
equations /,■ stands for the ligand activity in cell i and n is a 
measure of the amount of primary signaling source. Non- 
dimensional units are used with maximal signal production being 
1 for the primary signal without any loss of generality. 

We considered that each source can drive Notch activity 
production at a different rate such that it results into different 
values of the stationary saturated Notch activity. Parameter e 
accounts for the ratio of the stationary saturated Notch activity 
driven by cis-interactions over that one driven by the primary 
signaling source. A more complex biochemical reaction-based 
scheme indicates that t can be understood as the ratio between the 
signaling efficiency of the two sources, which depends on the 
signaling rate and the stability of each source (Methods). 
Therefore, we refer to t as relative signaling efficiency. t=0 
corresponds to the well-known cis-inhibition, in which cis- 
interactions titrate the receptor and drive no signaling. When 
f < 1 , cis-interactions drive signaling less efficiently than the 
primary signaling source. 

When the primary signaling source is driven by trans- 
interactions we have ii = r t (J i s ), where <(/,) stands for the weighted 
average of ligand activity in cells adjacent to cell i (Methods). 
Parameters r t c arising for each Notch signaling source are the 
trans and cis-interactions strengths, which are related to the trans 
and cis-binding and unbinding rates of the receptor and ligand 
complexes (see Methods for their definition in a more biochemical 
reaction-based approach). The trans and cis-interactions strengths 
parameters r l c set the threshold values of ligand for signal 
activation. 
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We set production of ligand to be inhibited by Notch signaling 
and took linear degradation for both signal and ligand activities 
(Methods). Taken together, the dynamics in cell i of the Notch 
signal Sj and the ligand /,- activities read: 



that cis-inhibition occurs for (Fig. 2D, Methods) 



dt 



-Pu+Pii- 



(3) 



dt 



\+bs'' 



(4) 



where time t is non-dimensional, v = vi/v s , being v/ and v s the 
degradation rates of the ligand and Notch signal respectively, b 
stands for the strength of ligand inhibition through Notch signaling 
and h represents an effective cooperativity of such an inhibition. 
When cis-interactions are not present (r c = 0), Eqs. 1-4 reduce to 
the model for Delta/Notch-mediated lateral inhibition dynamics 
early proposed by Collier et al. [13]. 

A switch from cis-activation to cis-inhibition 

We first evaluated whether the total Notch signaling increases 
(cis-activation) or decreases (cis-inhibition) when cis-interactions 
that drive Notch signaling are included. From Eqs. 1-3, it is 
obtained that cis-inhibition always occurs when cis-interactions do 
not drive a signal (f = 0) as expected (Fig. 2 A). In contrast, either 
cis-inhibition or cis-activation can arise when cis-interactions drive 
signaling on their own (Fig. 2B,G). 

We explored under which conditions each cis-signaling 
regulatory role (cis-inhibition versus cis-activation) arises and found 



1- 



(5) 



Otherwise, cis-interactions drive cis-activation (Fig. 2D). Accord- 
ing to the above relation, cis-inhibition requires cis-signaling to be 
less efficient than the primary signaling source {i.e. e< 1). Yet, less 
efficient cis-signaling does not ensure cis-inhibition. The above 
relation shows that the regulatory role of cis-signaling depends on 
the amount of the primary signaling source (j£j as well. When the 
primary signaling source comes from trans-interactions, the 
regulatory role of cis-signaling depends on the trans-interactions 
strength (r<). For e<l, the regulatory role switches from cis- 
activation to cis-inhibition as the primary source becomes more 
abundant (fi or r t increases) (Fig. 2B). 

The above result indicates that whether cis-signaling is acting as 
cis-activation or cis-inhibition does not depend on the strength of 
cis-interactions (r c ; Fig. SI). This is because the qualitative change 
(decrease or increase) in Notch signaling driven by the addition of 
ligand within a cell with primary Notch signal activity is 
independent of the amount of ligand being added. In contrast, 
the quantitative change depends on r c and on the amount of 
added ligand (Fig. SIC). 

Cis-signaling can drive two distinct effective circuit 
architectures 

The above analysis only took into account the dynamics of the 
signal when two sources of signaling (primary and cis-driven) are 
competing for the Notch receptor, Eqs. 1—3. We then asked which 
is the effect of this switch between cis-inhibition and cis-activation 
on the overall signal and ligand dynamics. To this end, we first 



B 





itr 3 ifr 2 icr 1 io° 




Figure 2. A switch between cis-activation and cis-inhibition. (A-C) Stationary Notch signal in a cell (.?,) versus the amount of ligand in that cell 
(/,) and the amount of primary signaling source (/() for (A) f = 0, (B) t = 0.1 and (C) e = 1. Red lines show the Notch signal dependence on /, in the 
absence of the primary signaling source (/i = 0) and for a primary signaling source with /( = 1. Decreasing curves indicate cis-inhibition and increasing 
curves show cis-activation. In B, cis-interactions drive cis-activation at low /< values, whereas they drive cis-inhibition at higher /(. (D) Parameter space 
showing where cis-activation (gray region) and cis-inhibition (white region) occurs, according to inequality 5. (E-F) Effective circuit architectures of 
the model when cis-interactions drive cis-activation (top) and cis-inhibition (bottom) for (E) isolated cells with a primary signaling source (straight 
arrow) and for (F) two adjacent cells that interact through trans-binding. Black arrows stand for activation, while red blunt arrows for inhibition. 
Parameter values: r c = 10 in panels A-C, h = 4 and i\ .= 1 in panel D. 
doi:1 0.1 371 /journal.pone.0095744.g002 
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evaluated the effective genetic circuits that arise when the 
dynamics of the ligand, Eq. 4, is included. By "effective" we 
mean that the circuit does not describe the individual interactions 
per se but their resulting regulatory role, which takes into account 
the context in which they occur. Therefore, we considered which 
circuit architecture arises when cis-interactions perform cis- 
activation and when they drive cis-inhibition. It can be readily 
seen that cis-interactions, coupled to a primary signal, give rise to 
two different effective genetic circuits (Fig. 2E). When cis- 
interactions drive cis-activation, a negative intracellular transcrip- 
tional feedback loop arises. In contrast, a positive feedback loop 
emerges when cis-interactions drive cis-inhibition. It is worth 
stressing that the amount of primary signal (jt) and the relative cis 
to primary signaling efficiencies (e) control which of the two 
effective architectures is acting by setting which is the regulatory 
role of cis-signaling. 

When the primary source is trans-interactions, these intracel- 
lular feedbacks are coupled to the intercellular mutual inhibition 
loop that is characteristic of lateral inhibition (Fig. 2F). Positive 
and negative feedback loops are well known to drive different 
dynamics, like bistability for the former and homeostasis for the 
latter (see for instance [55]). Accordingly, we can expect different 
roles of cis-signaling on patterning and cell fate choices, depending 
on whether it is in the cis-inhibition or cis-activation regime. 

Cis-activation inhibits fine-grained pattern formation, 
while optimal values of cis-inhibition promotes it 

From theoretical arguments it has been shown that the 
intercellular positive feedback mediated by trans-interactions is a 
sufficient mechanism for spontaneous pattern formation [13]. This 
feedback amplifies small differences in ligand and signal levels 
between cells and drives a mostly periodic lateral inhibition 
pattern composed of two cell types (Fig. 3A)[13]. This type of 
periodic pattern arises spontaneously for a large range of trans- 
interactions strengths (r t ) and above a minimal ligand inhibition 
strength (b) (Fig. S2 and Methods) [13]. The pattern solution exists 
and is stable in an even larger region of the parameter space, with 
a minimal yet lower ligand inhibition strength (Fig. S2 and 
Methods) [10,56]. 

From the effective circuit architectures (Fig. 2F), we may 
propose that cis-activation can inhibit pattern formation. Cis- 
activation drives a negative feedback loop within cells which could 
damper the amplification of differences between precursor cells 
driven by trans-interactions. In contrast, we may expect that cis- 
inhibition can promote patterning as it can enhance amplification 
of precursor differences by driving an additional positive feedback 
loop within cells (Fig. 2F). This latter expectation is in agreement 
with previous computational studies on cis-interactions that do not 
trigger a signal [25,29,33,35]. 

We first evaluated the case of efficient cis-signaling (t = 1) 
driving cis-activation whatever the strength and amount of trans- 
interactions. We explored extensively the parameter space to 
characterize where lateral inhibition patterning occurs (Methods). 
We chose a regime for which the lateral inhibition pattern can 
arise for a wide range of trans-interactions strengths (r r ) in the 
absence of cis-interactions (r c = 0; Fig. 3B). When cis-interactions 
are added (r c >0), pattern formation becomes forbidden (Fig. 3B). 
At very weak cis-interactions strengths, pattern formation is still 
possible albeit in a reduced range of trans-interactions strengths. 
We conclude that cis-activation inhibits patterning. 

Following the same procedures, we evaluated whether pattern 
formation is promoted by cis-inhibition. To this end, we analyzed 
the well-known case of cis-interactions which do not elicit signaling 
(fc = 0), driving always cis-inhibition. In this case, it is known that 




100 „ 200 300 
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Figure 3. Cis-activation inhibits patterning and cis-inhibition 
facilitates it. (A) Stationary lateral inhibition pattern formed in an array 
of irregular cells in the absence of cis-interactions (i ( =0). Grayscale is 
used to denote the ligand level (black for high ligand, /,■ = 1, and white 
for low ligand, /, = 0). (B-D) Regions of patterning for cis-interactions 
and trans-interactions strengths r c ,r, for (B) 6=1 (cis-activation), (C) 
e = 0 (cis-inhibition) and (D) e = 0.4. The black dot-dashed line in D 
divides the parameter space into the cis-activation region, on its left, 
and the cis-inhibition region, on its right. Blue regions show where the 
pattern grows spontaneously (LSA in Methods). Green dashed lines 
enclose the regions where lateral inhibition pattern solutions exist and 
are stable (Exact periodic solutions in Methods). In B, the patterning 
region is the one below the green dashed line. B and D show that 
patterning becomes forbidden as r c increases when cis-activation is 
acting. C and D show that patterning is enabled above a minimal cis- 
interactions strength r c when there is cis-inhibition. Other parameter 
values: h = 4 for all panels, r, = 1 and v= 1 for A, 6= 1000 for A-B and D, 
6 = 10 for C. 

doi:10.1371/journal.pone.0095744.g003 

cis-inhibition facilitates patterning by allowing it for graded trans- 
signaling and graded ligand inhibition (h = 1) [25] (Fig. S3A). We 
found that cis-inhibition facilitates spontaneous patterning as well 
for low ligand inhibition cooperativities (Fig. S3B). At higher 
cooperativities, the analysis showed that cis-inhibition promotes 
patterning too by reducing the minimal ligand inhibition strength 
(b) required for patterning (Fig. 3G). Cis-inhibition can have a 
detrimental effect as well (Fig. 3C). When the strength of cis- 
interactions is too high compared to trans-interactions strength 
[r c > > r t ), the coupling between cells becomes less relevant, 
impeding patterning. Yet, if the trans-interactions strength 
increases, patterning is enabled (Figs. 3C and S3A-C). 

We finally evaluated a more complex scenario in which cis- 
signaling switches from cis-activation to cis-inhibition when trans- 
interactions strength (r t ) increases (Fig. 3D). Our results show that 
in the cis-activation regime, pattern formation is inhibited, since 
an increase in cis-interactions impedes pattern formation. In 
contrast, in the cis-inhibition regime, pattern formation is 
facilitated since a minimal value of cis-interactions strength (r c ) 
enables spontaneous patterning. These results confirm the 
existence of distinct regulatory roles of cis-signaling (cis-activation 
versus cis-inhibition) as a function of the amount of trans- 
interactions as well as the different effect each of them has on 
lateral inhibition patterning. Taken together, the results suggest 
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that the effect of cis-signaling on patterning can be simplified to 
that of the regime in which it is acting. 

Cis-inhibition can modulate patterning and enhance 
multistability 

We next evaluated which patterning features arise in the cis- 
inhibition regime. Cis-inhibition can make cells worse receivers of 
inhibition [14,20,26,57]. This is confirmed in our model by 
evaluating the change in the threshold level of ligand activity 
required to drive ligand inhibition in an adjacent cell when cis- 
interactions are added in the receiving cell (Methods, Fig. S4). 
When cells become worse receivers of inhibition we can expect the 
ratio of high-ligand expressing cells to increase. Simulations results 
confirm cis-inhibition can increase the ratio of selected precursor 
cells (Fig. 4A for £ = 0 and Fig. S5 for e>0). The strength of cis- 
interactions (r c ) and the fraction between cis and trans-interactions 
strengths (r c /r,) become a control parameter for this ratio 
(Figs. 4B,C). When precursor cells exhibit large random initial 
variability between them in ligand and signal levels, the strength of 
cis-interactions can increase more gradually the ratio of high- 
ligand expressing cells (Figs. 4B-D). 

These results showed that cis-inhibition can enable a new 
regular salt-and-pepper pattern with 66% of cells highly expressing 



the ligand (Figs. 4A, S5). This pattern has the periodicity of the 
lateral inhibition pattern. However, the ratio of selected high- 
ligand expressing cells is complementary to it. Since our results 
show that cis-inhibition can modulate the threshold for lateral 
inhibition and thereby the ratio of selected precursors, we 
wondered whether it can enrich patterning and drive additional 
periodic patterns. By using a combined analytical-computational 
approach (Methods), we searched across the parameter space of 
cis-interactions (r c ) and trans-interactions strengths (r t ) whether 
and where different periodic patterns composed of two cell types 
were stable solutions of the dynamics. We chose to search for three 
different types of patterns that involve different numbers of 
selected precursors and spatial organizations (Fig. S6). One of 
them is the salt-and-pepper pattern of Fig. 4A with a 66% ratio of 
selected precursors (Fig. S6A). Another one is a salt-and-pepper 
pattern too but with a different periodicity and a 75% ratio of 
selected precursors (Fig. S6B). The third chosen pattern is stripped 
with 50% of cells expressing high-ligand levels (Fig. S6B). Our 
study showed that high enough cis-interactions strengths (r c ) 
enable the emergence of these patterns with high numbers of 
precursor cells (Figs. 5, S7, S8 and S9). In the absence of cis- 
interactions and for low cis-interaction strengths (r c —*0), the salt- 
and-pepper spatial organizations can persist but with much lower 




r c = 0 r c = 10 r c = 20 r c = 50 




r c = 0 r c = 10 r c = 20 r c = 50 



Figure 4. Cis-inhibiting interactions increase the ratio of high-ligand expressing cells. (A) Stationary patterns of ligand levels arising from 
precursor cells with small initial variability between them for different inhibiting cis-interactions strengths r c . Color code as in Fig. 3A. (B) Ratio of 
stationary high-ligand fated cells as a function of the cis-interactions strength r c when precursor cells show small (red triangles) and large (blue 
circles) initial variability between them. (C) Density plot representing the ratio of high-ligand cells in a tissue arising from precursor cells exhibiting 
large initial variability. Solid and dashed lines as defined in Fig. 3B-D. White vertical line is drawn for indicating the value c, = 3 along which 
simulations are performed in panels A, B and D. (D) Stationary patterns of ligand levels arising from precursor cells with large initial variability 
between them for different inhibiting cis-interactions strengths r c . In B-C panels, cells are considered high-ligand fated cells when its ligand level is 
over the threshold of 0.5. Parameter values: e = 0, /) = 4, b= 1000 and v = \ for all panels. Similar results are found for c = 0.4 (Fig. S5). In B, each point 
comes from the average of 10 numerical integrations of the dynamics on a lattice of 12 x 12 irregular cells starting at different initial conditions. In C, 
the results correspond to numerical integration of the dynamics performed over a lattice of 30 x 30 perfect hexagonal cells. 
doi:10.1371/journal.pone.0095744.g004 
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Figure 5. Cis-inhibiting interactions facilitate other periodic 
patterns to form. (A) Stationary stable stripped pattern of ligand 
levels that is a stable solution of the dynamics to small perturbations. 
Color code as in Fig. 3A. (B) Region (gray) where the pattern of stripes 
on a regular hexagonal array is a stable solution of the dynamics to 
small perturbations (Methods and Text SI) in the parameter space of cis 
and trans-interactions strengths r c and r t . Parameter values: /; = 4, e = 0 
and v=l for all panels and r, = 10 and r c = 100 for panel A. The stripped 
pattern appears also for e>0 in the cis-inhibition regime (data not 
shown). 

doi:10.1371/journal.pone.0095744.g005 

numbers of selected precursors (33% and 25%, Figs. S7, S8). In 
contrast, cis-inhibition with high enough cis-interactions strengths 
(r c ) enable the spatial organization of precursor cells within stripes 
(Figs. 5, S9). 

Additionally, we found that for high cis-interactions strengths 
(r c ) all these patterns are stable, i.e. there is multistability of pattern 
states (Figs. S7-S9). Hence, precursor cells could potentially 
become organized in any of them and should choose which 
pattern to form (Palau-Ortin et at, unpublished). 

Cis-inhibition allows localized patterning highly sensitive 
to the precursor state 

As shown in Fig. 4D, stable patterns without an obvious 
periodicity can arise too for high strengths of cis-inhibiting 
interactions (r c ). This occurs when precursor cells show large 
initial random variability between them, in agreement with [25]. 
We evaluated whether a high sensitivity to the initial state of 
precursor cells was causing the absence of periodicity. The results 
show that the finally formed stable pattern has strong similarity to 
the initial state of precursor cells (Fig. 6A-B). Simulations across 
the parameter space using precursor cells with large initial random 
variability between them confirmed that the pattern being formed 
is quite random as the state of precursors cells is for high cis- 
interactions strengths (Fig. S10). In contrast, large initial random 
variability between precursor cells dynamically evolves to more 
regular and periodic patterns of two cell fates for lower cis- 
interactions strengths (Fig. S10). 

Simulations results indicated that the patterns exhibiting high 
sensitivity to the initial random state of precursor cells keep 
spatially localized without spreading to the rest of the tissue 
(Fig. 6C). This absence of spreading is in sharp contrast with the 
dynamics of nucleating patterns driven only by trans-interactions. 
Nucleating patterns invade the rest of the tissue that is under 
lateral inhibition Notch dynamics through a traveling wave 
[58,59]. The localized patterns we find remain where they arise 
and do not spread despite the adjacent tissue is under lateral 
inhibition dynamics too. We searched across the parameter space 
of cis-interactions (r c ) and trans-interactions strengths (r t ) where 
pattern localization occurs (Fig. 6D). This pattern localization 
occurs in a specific region of the parameter space where cis- 
inhibiting interactions are dominant (r c /r t > 1). In addition, in this 



region there are many different stable pattern solutions and a 
homogeneous linearly stable state that impedes spontaneous 
patterning from small initial variability between precursor cells 
(Figs. 6D, S7-S9). These results suggest that strong cis-interactions 
performing cis-inhibition enrich patterning from precursor cells 
that show large initial variability between them. The arising 
patterns keep localized within the tissue and are reminiscent of the 
initial states of precursor cells. 

Cis-inhibition can drive cell-autonomous bistability 

We reasoned that cis-driven dynamics at the cell-autonomous 
level could be relevant for the phenomenon of localized 
patterning. Specifically, we wondered whether cis-inhibition could 
drive bistability of distinct ligand and signal level states in isolated 
cells. To evaluate it, we considered the role of cis-inhibition in the 
dynamics of single isolated cells that have a primary source of 
Notch signal (ji = fi 0 , being fi 0 a constant, see Text SI). This 
primary signaling source could be ligand-independent. 

Our results show that cis-inhibition can drive cell-autonomous 
bistability when a primary signaling source is present (Fig. 7A,B, 
Methods). This bistability drives similar cell types to the ones 
found in lateral inhibition patterning: cells are expected to be 
either on a high-ligand expression state or in a low-ligand 
expression state with opposite signaling state (Fig. 7A,B). Bistability 
of cell fates requires a minimal amount of cis-interactions (r c ) and 
of primary signaling source (ji; Fig. 7C). 

We evaluated the existence of this cell-autonomous bistability 
when the amount of primary signal corresponds to the signal that 
trans-interactions drive for the homogeneous state of equivalent 
cells (Methods). For this primary signal, bistability arises for high 
cis-interactions strengths (r c ) and encloses the region where 
patterns keep localized (Fig. 6D). This result suggests that cell- 
autonomous bistable dynamics arising from cis-inhibition may be 
relevant for the phenomenon of non-periodic pattern localization. 

Discussion 

Competition for signaling: a switch from cis-activation to 
cis-inhibition 

Several experimental evidences support the existence of 
signaling driven by cis-interactions [39-43], In this work we have 
theoretically characterized the effect of cis-signaling in different 
contexts. We found that a switch from cis-activation to cis- 
inhibition (or vice versa) arises. The switch can occur by 
quantitatively changing the signaling sources; either by changing 
the amount of the primary signaling source (e.g. trans-interactions), 
or by modulating the ratio between the signaling efficiencies of 
each source. As a result, phenotypes involving a reduction of 
Notch signaling when the ligand is increased within a cell (cis- 
inhibition) can be compatible with cis-signaling. 

The results show that cis-signaling can drive opposed capabil- 
ities to the patterning process, each arising on the different 
regulatory regimes of cis-activation and cis-inhibition. Cis-signal- 
ing acting as cis-activation creates a negative intracellular feedback 
loop that inhibits pattern formation. On the other side, cis- 
signaling acting as cis-inhibition creates a positive intracellular 
feedback loop facilitating patterning. This regime promotes 
patterning as cis-inhibition driven by null cis-signaling does 
[25,33,35]. 

The switch from cis-activation to cis-inhibition exemplifies a 
case of competition leading to a complex dynamical output: two 
signaling sources - the primary and the cis sources - with different 
efficiencies and competing for the same substrate - the Notch 
receptor - can result in a cis-productive signaling that drives cis- 
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Figure 6. Cis-inhibition allows pattern localization. (A-C) Initial (left) and stationary (right) patterns of ligand levels at high cis-interactions 
strengths r c in the cis-inhibition regime for different initial conditions: (A) all precursor cells have large initial random variability, (B) few precursor 
cells, distributed along a rectangle, have initial low ligand levels and (C) precursors within the top half of the tissue have initial high ligand levels and 
small variability, while precursors at the bottom half show large initial random variability in the level of ligand. In A-C, the final pattern strongly 
depends on the pattern formed by precursor cells. In (C) the pattern arises in a localized region (bottom half) and does not expand. (D) Region where 
localized patterns are found in a regular hexagonal array (gray) in the parameter space of cis and trans-interactions strengths r c and r,. Blue circles 
enclose the region for cell-autonomous bistability, where two states are linearly stable, according to simulation results (Methods). Solid and dashed 
lines as in Fig. 3B-D respectively. Parameter values: h = 4, 6 = 1000, 6 = 0 and v=l for all panels and r, = 3 and r c = 25 for (A-C). 
doi:10.1371/journal.pone.0095744.g006 



inhibition. Such effect is reminiscent of the behavior of full and 
partial pharmacological agonists, where a partial agonist can act as 
a competitive inhibitor of a full agonist [60] . It is also an example 
of complex regulation in which the interplay between different 
components changes the regulation performed by one of them 
[61]. _ 

This notion can be extrapolated to any competition between 
signaling sources that share the same receptor. When different 
ligands (canonical or not) bind the same type of Notch receptor but 
drive signaling with different efficiencies, one can expect 
competition between them and switches of regulatory roles [62] 
(Jelena et al., unpublished). This approach could also help to 
decipher the controversial roles of different non-canonical factors 
that bind to Notch and present both activatory and inhibitory 
effects on Notch signaling. For instance, this is the case of the 
Dlkl/2 non-canonical ligands [3,63] and the proteins MAGP1/2 
[45,64]. Different competition events have already been shown to 
drive significant regulatory effects in other signaling pathways (see 
for instance [65-70]). 

In addition, our approach provides a definition for the ratio of 
signaling efficiencies coming from different sources that share the 
same receptor. This ratio is defined by the signaling rates and by 
the stability of the signaling sources. This definition is relevant to 
determine the regulatory role each signaling source drives on the 
overall signaling. 

Cis-inhibition as a modulator of the ratio of selected 
precursor cells. From experimental grounds it has been 
already pointed out that cis-inhibition can drive cells to become 
worse receivers [14,20,26,57]. This effect can be obtained from 
our model too and enables cis-interactions to increase the number 
of high-ligand expressing cells. Cis-inhibition is expected also to 
drive cells to become worse signal senders by sequestering the 
ligand. Albeit this aspect is not considered in our simplified model, 
this effect should increase the ratio of cells reaching the high-ligand 
fate too. We have checked that a more complex model involving 
ligand sequestration similar to Sprinzak et al. (2010) (Methods and 



Text SI) also shows an increase in the ratio of high-ligand 
expressing cells. Consistent with our simplified model, competition 
between signaling sources in the complex model also yields 
switches of the regulatory role of cis-interactions (Fig. S 1 1 , 
Methods). The results confirm the increase in the ratio of high- 
ligand fated cells with the strength of cis-interactions in the cis- 
inhibition regime (Fig. SI 2). This complex model indicates that 
this increase occurs too when there is no cooperativity in the 
inhibition of the ligand {i.e. for h=\ in Eq. S4a of Text SI, Fig. 
SI 2). 

It has been reported that Lunatic Fringe knockdown produces 
an increase in the neurogenesis ratio in the hindbrain of zebrafish 
embryos [71]. This could be an example of the cell- type ratio 
modulation we find in our simulations, since Fringe potentiates 
Delta-Notch trans-interactions [72] and could inhibit the cell- 
autonomous association of Delta and Notch [20] . 

The strength of cis-inhibiting interactions is not the only 
potential modulator of the ratio of high-ligand expressing cells. 
Specifically, different theoretical approaches that do not take cis- 
interactions into account have reported other components that can 
modulate this ratio through changes in the level of the threshold 
that drives inhibition of the ligand [73,74]. Also, it has been shown 
that higher rates of Delta production can drive a graded increase 
of high-Delta cells, which was validated experimentally [75]. 

Noteworthy, Notch signaling dynamics in vivo can select a 
different ratio of high-ligand cells in different contexts [76-81]. 
Our results suggest that cis-inhibition could underlie the selection 
processes that involve high ratios of selected cells. 

Cis-inhibition potentiates multistability and enriches 
patterning. Multistability can enable the change of fate of cells 
and it has been widely evaluated in the context of single cells, 
specially in bacteria and stem cells. A recurrent circuit topology 
that shows multistability and participates in stem cell renewal and 
differentiation is the toggle switch with auto-activation [82-84]. 
Auto-activation in this dynamics facilitates multistability [82]. 
Herein we show that lateral inhibition Notch dynamics with cis- 
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Figure 7. Cis-inhibition with a primary Notch signaling source 
creates cell-autonomous bistability. (A) Stationary ligand level as a 
function of the cis-interactions strength r c for ^ = ^ 0 = 10 and t = 0. 
Solid lines denote linearly stable solutions, dashed lines indicate linearly 
unstable solutions. Black dots refer to the stationary ligand levels for 
JV = 100. (8) Nullclines diagram showing the three possible solutions at 
r c = 100. The blue and red lines represent the nullclines. The continuous 
black line is a separatrix, which divides the parameter space into two 
basins of attraction of the two stable solutions. Percentages indicate the 
fraction of cells reaching the corresponding stable state computed from 
4xl0 4 cells with initial random uniform levels of ligand. (C) Phase 
diagram showing the cell-autonomous bistability region zone where 
two states are linearly stable. The gray area is the theoretically 
computed region, and the blue circles correspond to simulation results 
(Methods). Parameter values: h=4, 6=1000, t = 0 and v = l for all 
panels. These results can also be obtained for f >0 in the cis-inhibition 
regime (data not shown). 
doi:1 0.1 371 /journal.pone.0095744.g007 

inhibition can be described with an effective topology that 
corresponds to a toggle switch with auto-activation (Fig. 2F). In 
this effective topology, cis-inhibition drives auto-activation and 
facilitates multistability, in agreement with the effect of auto- 
activation in cell-autonomous toggle switches. This is reminiscent 
to the reported multistability in the ommatidia formation in 
Dro.sophila eye, which has been proposed to be driven by an auto- 
activatory feedback loop due to Atonal [85]. 

At the cell-autonomous level, we find that cis-signaling can drive 
bistability of ligand and signaling states when a basal cell- 
autonomous activity of Notch is present and cis-inhibition is 
acting. In this case, the effective circuit topology corresponds to a 
positive feedback loop that involves a mixed-feedback loop 



[86,87]. The bistability regime is confirmed and becomes more 
prominent in the more complex model (Fig. SI 3, Methods). The 
prediction of cell-autonomous bistability due to cis-inhibition 
could shed light to new functions of Notch in single cells. Recently, 
cell-autonomous bistability in Notch has been identified in the 
context of colon cancer stem cells [88] . In particular, it has been 
shown that the sequestering of mRNA Notch 1 by the tumor 
suppressor microRNA miR-34a drives cells with Notch signal 
bimodality [88]. 

Simulation results show that cis-inhibition enables the spatial 
localization of patterns, which do not propagate spontaneously on 
the entire tissue. This could correspond to a wave-pinning 
phenomenon [58]. Typically, wave-pinning arises in discrete 
dynamical systems when the coupling between the discrete units is 
below a critical strength [89]. In our scenario, the critical coupling 
would be related to cis versus trans-interactions strengths (r c /r t 
ratios). Moreover, the most disordered patterns, with high-ligand 
cells adjacent to each other [25], appear in the region where 
localized patterning occurs. In such regions of the parameter 
space, the final pattern strongly depends on the initial precursor 
state. This dependence on the precursor state is reminiscent of the 
directionality provided by cis-inhibition in the differentiation of 
R1/R6/R7 precursor photoreceptor cells in the Drosophila eye 
[26]. Together, these results suggest that cis-inhibition can enrich 
patterning by enabling additional modulations of cell fate 
decisions. 

Methods 

Model formulation within an irregular cellular array 
layout 

The simple model phenomenologically includes the competition 
between cis and the primary signaling sources for the Notch 
receptor and the inhibition of Notch signaling on the ligand. It is 
based on the approach introduced by Collier et al. [13] for lateral 
inhibition dynamics through trans-interactions. We included 
competition such that it is in agreement with a more biochemical 
reaction-based approach (see Complex model below). The 
simplicity of the simple model strongly facilitates the vast 
exploration of different patterning regions in the parameter space. 

The weighted average of non-dimensional ligand concentration, 
</,->, appearing in Eq. 1 through /( = r,</,> due to trans- 
interactions, describes the interactions between adjacent cells on 
a two-dimensional irregular array of cells: 



</<•> = 



(6) 



being e,y the length of the cell membrane 



witha,y = e,://E /e „„ 
edge shared by adjacent cells i and j, and the summation involves 
all cells adjacent to cell i (jenn(i)) [90]. We constructed an 
irregular two-dimensional array of cells with periodic boundary 
conditions as in [10] with irregularity parameter y= 1.33 (see Fig. 
S3 in [10]). First, we generated an irregular distribution of points 
on a plane starting from a perfect triangular lattice and considered 
periodic boundary conditions by surrounding the array of points 
with equivalent arrays. Second, a Voronoi tessellation was created 
around these points using Mathematica's Computational Geom- 
etry Package (Wolfram Research, Inc. (2008), Mathematica, 
Version 7.0, Champaign, IL, USA). 
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Complex model 

The Complex model takes into account the dynamics of the 
Notch receptor and of the complexes formed by receptors and 
ligands (see Text S 1 for all reactions, model equations and further 
details). It includes receptor and ligand inactivation through 
proteolytic cleavage [25,35] and it does not make assumptions 
regarding the capability of sending and receiving signals (e.g. when 
cis-interactions are acting, signal sending cells are not necessarily 
refractory to receive inhibitory signals from its neighbors). In this 
Complex model, the dynamics of the free receptor (TV,), trans and 
cis-formed receptor-ligand complexes (C,- and B, respectively), and 
Notch signal «S, in cell i when the primary signaling source is due 
to trans-interactions read: 

A AT. 

-r± = -k b N i <.L i }-k' b N i L i +k- b C i +k'_ b B i + IS n - Vn N i (7) 



and e= /- l f-J + Vc . Ln is a characteristic dimensional concentration 

of ligand (i.e. Lj = Ldlj). The first term on the right-hand side 
corresponds to the stationary primary signaling driven by trans- 
interactions whereas the second term is the stationary signaling 
driven by cis-interactions. 

We define the efficiency of each source as the ratio of success to 
signal of the receptor-ligand complexes. This efficiency corre- 
sponds to k s /(k s + v c ) for the primary signaling source and to 
k' s /{k' s + v b ) for cis-interactions. Therefore, £ parameter 

( e= FF+vp ' s me re l a tive efficiency of the cis-driven source 
compared to that of the primary signaling source. From Eq. 1 1 it 
can be seen that e corresponds as well to the ratio of maximal 
saturated stationary Notch activity driven by cis-interactions over 
that one driven by the primary signaling. 

The equations of the Complex model for single isolated cells 
with ligand-independent and cell-autonomous primary signaling 
sources are detailed in Text SI. 



Evaluation of the regulatory role of cis-interactions 

We defined the regulatory role of cis-interactions (cis-inhibition 
— j— =k b Ni(Li) — k- b Ci — k s Ci — v c Ci (8) or cis-activation) through the (negative or positive, respectively) 

change in Notch signal dynamics within a cell when its ligand 
content increases, R Sl ,l,'. 



dS, 

~d7 



=k' b N i L i -k'_ b B i -k> s B i -v b B i 



(9) 



d7 



= k s Q + k'B i -v s S i 



(10) 



where the variables are in dimensional units and Lj is the free 
ligand in cell i, whose dynamics are detailed in Text SI. is 
given by Eq. 6 applied on species L. % is the dimensional time. The 
binding and unbinding dynamics of Notch receptors with its ligand 
in trans and in cis have rates k b , k_ b and k' h , k!_ b , respectively. 
Notch production (fi n ) and degradation (v„) are also taken into 
account. Trans and cis complexes, C; and 2?,-, have degradation 
rates v c and v b , respectively. The model does not detail the overall 
mechanism by which trans-interactions drive signal activity. 
Instead, it assigns a rate k s to the proteolytic cleavage of the trans 
complex and the ultimate release of Notch signal S;. The model 
also considers the case in which cis-interactions drive Notch 
signaling. We implemented it by taking into account the argued 
mechanism for cis-signaling [40], so that the release of Notch 
intracellular domain would also occur for cis complexes (Bj). We 
set this step to occur at rate k' s . By taking k' s = 0, the above 
equations account for the usual scenario of cis-interactions that 
sequester the receptor and drive no signaling. 

Notice that the stationary solution of Eqs. 1-3 is the same 
function of and (Ji) as the stationary solution of Eqs. 7-10, 
which reads (diV,/dT = 0, dC,/dr = 0, dS,/dT = 0, dS,/dT = 0): 



C st_r.st r t(Ji) i est £r Ji 



with Sjf = 



vj(ij+v c )' 



vji k_f ! +ks + v c 



Lo, r c 



R s i = 



j^dsj_ _ r c (e + fi(e-l)) 
dk At ~ (l+n+rjj) ' 



(12) 



where the result of the derivative for the model described by Eqs. 
1-3 is indicated. Cis-inhibition is defined as a decrease in Notch 
signaling when the ligand content increases within the same cell, 
i.e. RsiJj < 0, whereas cis-activation corresponds to an increase in 
Notch signaling, i.e. R Si ,Ij > 0. Based on the above expression for 
Rs-,1-, cis-interactions drive cis-inhibition when e + fi(e—l)<0. 
Therefore, the condition for cis-inhibition to happen can be re- 
written as inequality 5: 



£< 



l+H l+r ( </i>' 



(13) 



where H = r t (li) for trans-interactions has been introduced in the 
last right-hand side term. This inequality states that cis-inhibition 
occurs when the maximal (saturated) signaling driven by cis- 
interactions is lower than the signaling driven by the primary 
source when acting alone (Fig. SI). Notice that this criterion for 
the cis-regulatory role is independent of r c (Fig. SIC). 

Eq. 13 gives the regulatory role of cis-interactions in the 

Complex model at the steady state as well (Rs,,L, = 757) with 



— k,+Ve ->-"A r -—^r k ^iJ' +v Lq. Accordingly, for no cis- 



signaling (k' s = 0) cis-interactions perform cis-inhibition, whatever 
the context and additional parameter values. For k! s jk s > v b /v c (i.e. 
e > 1), cis-interactions perform always cis-activation. For 
0<k' s /k s <v b /v c (i.e. 0<e<l) a switch from cis-activation to 
cis-inhibition can occur as the amount of trans-interactions 
increase. Noteworthy, when the receptor-ligand complex formed 
by cis-interactions is more unstable than the complex formed by 
trans-interactions (v b /v c >l), it can drive cis-inhibition even if it 
signals faster than the trans complex (1 <k'Jk s <v b /v c ). 

From the model equations, it can be readily seen that the switch 
of regulatory role can only take place in the non-linear regime of 
the signaling function (Eqs. 1, 2 and 11). Notice that this regime 
does not require saturation of the Notch receptors. When the 
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primary source is acting in the linear regime (i.e. fi < < 1), the 
addition of ligand within the cell does not reduce the primary 
signaling since there is no competition for the Notch receptor. As a 
result, cis-signaling always drives cis-activation in this linear 
regime. 

The inequality arising for isolated cells with cis-signaling and a 
primary signaling source is detailed in Text SI (Fig. SI IB). 



values (dsi/dt = 0, d/,/d? = 0) of these two cell types are: 



1 



r t <.I A } + er c l A 



(16) 



Linear stability analysis (LSA) 

Linear stability analysis [13,56] has been applied to Eqs. 1-4 to 
determine in which regions of the parameter space spontaneous 
patterning occurs. We defined spontaneous patterning as the 
process that drives pattern formation from a linear instability of 
the homogeneous initial state through small non-homogeneous 
perturbations (i.e. when small initial variability between precursor 
cells becomes amplified) [91]. LSA enabled us to make analytic 
predictions of how the pattern formation capabilities of the system 
would be changed by cis-interactions. LSA indicated that cis- 
interactions do not change the fastest growing mode, and hence 
the periodicity of the pattern is expected to be the same as in the 
absence of cis-interactions for spontaneous patterning (Text SI). 
LSA indicated that spontaneous patterning in a regular hexagonal 
array of cells with periodic boundary conditions would happen 
when (see Text SI for details and Fig. SI 4) 



l<^(W-sf/2), 



being 



IB^dh 
v dsi dt 



Vo 



* dh dt 



Vo 



d dsi 
8l k dt 



(14) 



(15) 



•Vo 



where 4 is the ligand level in a neighboring cell to cell i, so and /o 
are the homogeneous steady states for the multicellular system (i.e. 
the solutions of d/,/dr = 0, ds,/dt = 0 for </,> = /,) and tu is the 
number of nearest neighbors in a hexagonal cellular array (to = 6). 
3§ measures the strength of ligand repression at the homogeneous 
stationary state and verifies £$ < 0. sf measures the strength of 
trans-activation and verifies ,$/>0. Notice that W is R,. /. 
computed at the homogeneous steady state. Therefore, when cis- 
inhibition (R s .j. < 0) is acting at such homogeneous state then we 
have W<0. W measures the strength of cis-inhibition (when W<0) 
and of cis-activation (when g 7 >0). 

According to inequality 14, cis-inhibition (e?<0) facilitates 
patterning (Fig. SI 5). In contrast, cis-activation (^?>0) inhibits 
patterning (Fig. SI 5). Inequality 14 was also used to evaluate 
where spontaneous patterning can emerge in the r t -r c parameter 
space. These results are depicted by solid lines in Figs. 3, 4 and 6. 
The use of the simple model of Eqs. 1-4 enabled a vast exploration 
across the parameter space. We checked several of the regions 
obtained by LSA with numerical simulations (Text SI, and Fig. 
S16 as example). 

Exact periodic solutions 

We evaluated the r t -r c parameter space regions where the 
lateral inhibition pattern is a stable solution of the dynamics 
defined by Eqs. 1-4. We also evaluated whether other periodic 
patterns are stable solutions of these dynamics. This analysis was 
strongly facilitated by the use of the simple model. To this end, we 
extended a method we previously introduced [56] to our system 
and to new periodic patterns. We considered periodic patterns 
composed of only two different cell types: cell type A and cell type 
B. According to dynamics given by Eqs. 1-4, the stationary state 



1 

1+bA' 



sb-- 



r t <lis> + £r c l A 
\+r,{l B s ) + r c l B 



(17) 



Based on the periodicity of the pattern, we imposed which is the 
neighborhood of cell types each cell type interacts with: 

(J/)=m A l A +{\-m A )h, (lBy=in B lB + (l-mB)l A , (18) 

where m AtB is the ratio of ^4,i?-like cells neighboring to the A,B- 
cell type. For the common lateral inhibition pattern we have 
m A = 0 and mg = l/2. The r t —r c parameter boundary regions 
enclosing the region where this lateral inhibition pattern is a stable 
solution of the dynamics are depicted with dashed lines in Figs. 3, 
4 and 6. The m A ,B values for other periodic patterns are as follows. 
For an additional salt-and-pepper pattern (Fig. S6B) we have 
m A =Q and rag = 2/3. For the stripped pattern m A = l/3 and 
wig = l/3 (Fig. S6B). 

Solutions of Eqs. 16-18 were found using NSolve from 
Mathematica and also through custom made programs using the 
bisection method. Stability of solutions was evaluated computa- 
tionally by numerical integration of the dynamics with patterned 
initial conditions (see below and Text SI). Together, these results 
show that the parameter space region where the pattern with the 
periodicity of the lateral inhibition pattern is stable is the largest 
one and contains the regions where the other patterns are stable 
(Figs. S7-S9). 

Threshold for lateral inhibition 

We evaluated how cis-interactions within a cell (i) change its 
capacity to receive the inhibition from adjacent cells. We termed /* 
the threshold for lateral inhibition and defined it as the ligand 
activity in adjacent cells that drives the inhibition of the ligand in 
cell i. This inhibition of the ligand was defined as having a 
production rate of ligand activity x times (0<x< 1) the maximum 
production rate, which is 1 for Eq. 4. Therefore, according to Eq. 
4, a production rate x of ligand occurs for s* =((1 — x)/(xb)f i ^ , \ 
We then computed which ligand activity in adjacent cells /* is 
required to drive a stationary signaling s* within cell i if this cell 
has a ligand activity /; which drives cis-interactions (i.e. we isolated 
/* from s* = i+rj,"'' 1 1. according to Eqs. 1-3). Taken together we 
obtain that the threshold for lateral inhibition is: 



/* = /* 



1+^(1-^ 



(19) 



where It 



is the threshold for lateral inhibition in the 



absence of cis-interactions. The above equation indicates that cis- 
interactions increase this threshold (/* > l* rans ) making cells worse 
receivers of inhibition when cis-inhibition is taking place (Fig. 4A- 
C). 
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Cis-driven cell-autonomous bistability 

We evaluated whether the positive feedback generated by cis- 
inhibition with a competing primary signaling source in isolated 
cells (Fig. 2E) is sufficient to drive cell-autonomous bistability. To 
this end, we computed the steady state solutions (d.s,/d? = 0, 
d/,/d? = 0) for Eqs. 1-4 when H = fi 0 . Linear stability of these 
solutions was given by 

\<£3W, (20) 

where 33 and are defined as in the LSA section in Methods but 
are evaluated at the steady state solutions. Analysis of the nullclines 
(d.?,/d? = 0, d/,/d? = 0) indicates that cis-inhibition is required to 
have bistability (both nullclines need to be decreasing functions). 
We extensively explored in the \i$-t c parameter space where 
bistability occurred with LSA, and corroborated it with numerical 
simulations (Fig. 7C). In Fig. 6D we used fi 0 = r t lo being /o the 
homogeneous steady state for the multicellular system. All these 
simulations were performed for 400 cells with random initial 
conditions. The region of bistability in the phase space was 
delimited similarly than the LSA regions (Text SI). The stability of 
the bifurcation branches in Fig. 7A was checked with simulations 
of 900 cells with small fluctuations around the different branches. 

Cell-autonomous bistability in the Complex model was 
analyzed through nullclines analysis and numerical integration of 
the dynamics. In this case, bistability arises too in the absence of 
cooperativity in ligand inhibition (h = l). 

Fixed points were computed with Mathematica with the NSolve 
function and with custom made programs with a bisection 
method. 

Numerical integration of the dynamics 

We integrated Eqs. 1-4 with custom made programs using a 
Runge-Kutta fourth-order algorithm [92] with a time step of 
d? = 0. 1 . The Complex model for the multicellular and single cell 
system was integrated with Mathematica by using the NDSolve 
function. 

To evaluate which patterns were formed and their stability we 
used three types of initial conditions. (1) Precursor cells (i.e. cells at 
their initial condition) show small variability between them: for 
each molecular species x, x,(t = 0) = (7,(1 + d(Uf —0.5)) being u, 
the homogeneous steady state for the specie x, Uf a uniform 
random number between 0 and 1 and (5 = 0.1. (2) Precursor cells 
show large random variability between them: x,-(f = 0)= Uf, 
unless otherwise stated. (3) Precursor cells already form a regular 
pattern, with small variability: X/(t = 0) = £7,(1 + S(Uf — 0.5)) being 
c, the steady state pattern solution and 8 = 0.1. We performed 
numerical integration of the dynamics for different parameter 
values defined on a logarithmic mesh across the r t -r c parameter 
space. Simulations were stopped when the steady state was 
reached. 

Supporting Information 

Figure SI The regulatory role of cis-interactions when 
acting together with trans-interactions. (A C) Signal 
production rate in cell i due to (left) trans P\, and (middle) cis- 
interactions Pn, and (right) total signal production rate Pi; + P%i as 
a function of the ligand level /, within the cell. (Left) P\i decreases 
with for all parameter values because cis-interactions drive 
competition for the Notch receptor. </,-> value is set at the 
homogeneous fixed point. (Right) P% always increases with /, since 
we impose that cis-signaling on its own activates Notch signaling. 
(Right) Pn + Pii can be either a decreasing or an increasing 



function of/,-. This indicates the regulatory role of cis-interactions. 
Cis-inhibition occurs when Pu + Px decreases with /;, whereas cis- 
activation is acting when Pu + Pli increases with /,. (A) Signal 
production rates for different ratios between cis/trans signaling 
efficiencies c. £ = 0.1 (gray solid line), £ = 0.5 (black solid line) and 
£ = 0.9 (dashed line). The regulatory role changes from cis- 
inhibition to cis-activation as £ increases. (B) Signal production 
rates for different trans-interactions strengths r t : r, = 50 (gray solid 
line), r, = 100 (black solid line) and r, = 1000 (dashed line). The 
regulatory role changes from cis-activation to cis-inhibition as r, 
increases. (C) Signal production rates for different cis-interactions 
strengths r c : r c = 10 (dashed line), r c = 100 (black solid line) and 
r c = 250 (gray solid line). The regulatory role does not change with 
r c . This is because neither Pn at /, = 0 
(Plitt' = 0) = r,</,->/(l +r,</,») nor the saturated value of Pjx, 
P^f = e, depend on r c . Parameter values are h = 2, 6=1000, 
r, = 250, r c = 500 if not indicated otherwise, (B) £ = 0.5 and (C) 
£ = 0.9. 
(TIFF) 

Figure S2 Results in the absence of cis-interactions 

(r c = 0). (A) Scheme of interactions as in [13] of two cells that 
inhibit each other through Notch-mediated lateral inhibition. 
Black (blunt red) arrows denote activation (inhibition). Notice the 
positive intercellular feedback loop. (B-C) Phase diagrams in the 
parameter space of ligand inhibition strength b and trans- 
interactions strength r t for (B) high (h = 4) and (C) low (h = 2) 
cooperativity in ligand inhibition. The blue region in (B) is where 
the homogeneous state is linearly unstable. This is the region of 
spontaneous patterning, where the lateral inhibition pattern can 
arise from the amplification of small differences between precursor 
cells, as described in [13]. The region above the dashed line is 
where the pattern solution (with the periodicity shown in Fig. 3A) 
is an exact stable solution of the dynamics [56]. Above the dashed 
line and below the blue line in panel B both the homogeneous 
state and the lateral inhibition pattern are stable solutions of the 
dynamics (i.e. it is a bistable region). The continuous and dashed 
lines in (B) have been shown in [10]. Spontaneous patterning does 
not occur at low cooperativities (h = 2). 
(TIFF) 

Figure S3 Phase diagrams in the r t —r c parameter space 
for cis-inhibition with different cooperativities in ligand 
inhibition. (A— C) Phase diagrams in the r t -r c parameter space 
for (A) no cooperativity (h= 1), (B) low (h = 2) and (C) high (A = 4) 
cooperativity. (A) In the absence of cooperativity (h = 1), a minimal 
amount of cis-interactions is required to create a pattern for any b 
value, being consistent with Sprinzak et al. (2010) [25]. (B) At low 
(h = 2) cooperativity, cis-interactions enable spontaneous pattern- 
ing. (C) At high cooperativity (h = 4) cis-interactions can promote 
the bistable regions. In all panels, b= 1000 and £ = 0. Color codes 
and line types as in Fig. S2. 
(TIFF) 

Figure S4 Cis-interactions in the cis-inhibition regime 
make cells worse receivers of inhibition. (A) Threshold for 
lateral inhibition /* (Eq. 19 with x = 0.01 and /, = 0.5) as a function 
of the trans-interactions strength r t for £ = 0. Results for different 
cis-interactions strengths are depicted: r c = 0 (solid line), r c =W 
(gray dashed line) and r c = l00 (dotted-dashed gray line). The 
vertical line is a guide to the eye for a particular trans-interactions 
strength value, to better appreciate the rise of /* due to cis- 
interactions strength. (B-C) Contour lines for different /* values 
are depicted across the r t -r c parameter space for (B) e = 0 and (C) 
£ = 0.4. Lines are depicted for /*=0.1 (long-dashed), 0.5, 0.7, 1.0 
(short-dashed). As a guide to the eye, the spontaneous pattern 
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formation regions (enclosed by blue lines) and the regions where 
the pattern is a stable solution of the dynamics (enclosed by green 
dashed lines) are depicted. Other parameter values are as in Fig. 
S3C. 
(TIFF) 

Figure S5 Cis-inhibiting interactions increase the ratio 
of high-ligand expressing cells at £ = 0.4. Simulation results 
showing patterns of ligand levels from precursors with large initial 
variability between them for different cis-interactions strengths (r c ). 
Grayscale is used to denote the ligand level (black for the highest 
ligand activity, /, = 0.035, and white for no ligand activity, // = 0). 
Other parameter values are r, = 1000, b= 1000, h = 4 and v= 1. 
(TIFF) 

Figure S6 Representation of periodic patterns com- 
posed of two cell types on a regular hexagonal array. (A) 

Salt-and-pepper patterns with the periodicity of the fastest growing 
mode (Sl = — 1/2, see LSA in Text S 1). The P and I patterns have 
the same periodicity but 33% and 66% of cells, respectively, are 
high-ligand expressing cells (black). (B) Patterns with the 
periodicities of the secondary fastest growing modes (fi = — 1/3, 
see LSA in Text SI). P2 and 12 are salt-and-pepper patterns too 
with 25% and 75% of high-ligand expressing cells respectively. 
The pattern of stripes (S) has 50% of cells with high ligand levels. 
On the right of each row of patterns, two groups of 6 cells 
neighboring a central cell illustrate how many neighboring cells 
are like the central one and how many are different. Each group 
has a different cell type on the center (cell type A in violet and cell 
type B in green). Notice that cell types A and B are defined by the 
rriA and mg values (Methods) and not by their ligand level. These 
illustrations facilitate the computation of itia and ms values of Eqs. 
16-18 for each pattern. 
(TIFF) 

Figure S7 Cis-inhibiting interactions enable the salt- 
and-pepper pattern with 66 % of cells highly expressing 
the ligand. (A) Phase diagram showing where patterns (green) P 
and (red) I (as defined in Fig. S6) are each a stable solution of the 
dynamics to small perturbations. Green dashed and blue lines as in 
Fig. S3C. (B) Bifurcation diagrams for each cell type, A and B, for 
r t = 3. The periodic solutions are shown in blue. The homoge- 
neous solution is shown in gray. Solid (dashed) lines correspond to 
linearly stable (unstable) states. At low cis-interactions strengths, 
the stable branches correspond to P (light blue) and at higher cis- 
interactions strengths to I (dark blue). Note that there is a large 
parameter region in which both patterns are stable. Solutions for 
patterns were found by solving Eqs. 16-18 with m J 4=0 and 
wig = 1/2. Stability of solutions was evaluated through numerical 
simulations of the dynamics (Text SI). Parameter values: £ = 0, 
h = 4, 6=1000 and v=l. 
(TIFF) 

Figure S8 Cis-inhibiting interactions enable the salt- 
and-pepper pattern with 75% of cells highly expressing 
the ligand. (A) Phase diagram showing where patterns (green) P2 
and (red) 12 (as defined in Fig. S6) are each a stable solution of the 
dynamics to small perturbations. Green dashed and blue lines as in 
Fig. S3C. (B) Bifurcation diagrams for each cell type, A and B, for 
r, = 3. The periodic solutions are shown in blue. The homoge- 
neous solution is shown in gray. Solid (dashed) lines correspond to 
linearly stable (unstable) states. At low cis-interactions strengths, 
the stable branches correspond to P2 (light blue) and at higher cis- 
interactions strengths to 12 (dark blue). Note that there is a large 
parameter region in which both patterns are stable. Solutions for 
patterns were found by solving Eqs. 16-18 with niA=0 and 



wig = 2/3. Stability of solutions was evaluated through numerical 
simulations of the dynamics (Text SI). Parameter values: £ = 0, 
h = 4, 6=1000 and v=l. 
(TIFF) 

Figure S9 Cis-inhibiting interactions enable the 
stripped pattern with 50% of cells highly expressing 
the ligand. (A) Phase diagram showing where pattern (green) S 
(as defined in Fig. S6) is a stable solution of the dynamics to small 
perturbations. Green dashed and blue lines as in Fig. S3C. (B) 
Bifurcation diagrams for each cell type, A and B, for r t = 3. The 
periodic patterned solution is shown in blue. The homogeneous 
solution is shown in gray. Solid (dashed) lines correspond to 
linearly stable (unstable) states. Since type A and type B cells are 
equivalent for this pattern, there is bistability of the stripes 
solution. In all the parameter region where the stripped pattern is 
stable there are several patterns (P, I, P2 or P2) that are stable too 
(Figs. S7, S8). Solutions for patterns were found by solving Eqs. 16- 
18 with m^ = l/3 and nig =1/3. Stability of solutions was 
evaluated through numerical simulations of the dynamics (Text 
SI). Parameter values: £ = 0, h = 4, b= 1000 and v= 1. 
(TIFF) 

Figure S10 Sensitivity to initial conditions occurs at 
high cis-interactions strengths in the cis-inhibition 
regime. (A— C) Stationary patterns of ligand level emerging from 
precursors with large initial random variability between them. (D- 
F) Structure function (Text SI) of the patterns in A-C respectively, 
without the homogeneous mode (SI = 1). (A,D) The disordered 
pattern that emerges at high cis-interactions strengths (r c = 1000), 
where the homogeneous solution is linearly stable. (B,E) Regular 
pattern that emerges at lower cis-interactions strengths (r c = 200), 
where the homogeneous solution is linearly stable. In these panels 
we set r t = 50. (C,F) Salt-and-pepper periodic pattern that emerges 
within the spontaneous pattern formation region, for very low cis- 
interactions strengths (r c = 0A and r ; = 100). Other parameter 
values are £ = 0, h = 4, b= 1000 and v= 1. 
(TIFF) 

Figure Sll The switch between cis-activation and cis- 
inhibition regulatory roles also occurs in the Complex 
model. Stationary Notch signal in cell t, Eq. 1 1 , versus the amount 
of free ligand in the cell, L, , and the primary signaling source for 
(A-C) the multicellular system (K,<L,> with K t = kt,/k_i,) and (D- 
F) the single cell system (Kq = jU 0 /j«-o) for ( A > D ) nul1 iK = % ( B > E ) 
slow (k' s = 0Ak s in B, and ^ = 0.05^ in E) and (C,F) fast k' s = k s 
cis-signaling. The value of £ is (A,D) £ = 0, (B) £ = 0.4, (E) £ = 0.24 
and (C,F) £= 1. Red lines show the dependence of Sj on L, when 
there is no primary source and when it is maximal on the plot. An 
increasing function denotes cis-activation, while a decreasing 
function corresponds to cis-inhibition. A,D (£ = 0) show cis- 
inhibition; B,E (£ = 0.4 and £ = 0.24) show a switch from cis- 
activation to cis-inhibition as the primary source increases; D,F 
(e = 1) show cis-activation. Other parameter values: /?„ = 1 au 
hr- 1 , k' h = \0 au" 1 hr- 1 , k'_ h = 0A hr" 1 , v c = v A = 0.1 hr" 1 , 
v„ = l hr -1 and v, = k s = 0.5 hr -1 for all panels; = 0.1 hr -1 
for (A-C) and ji_ 0 = 0.1 hr -1 for (D-F). hr refers to hours and au 
refer to arbitrary concentration units. 
(TIFF) 

Figure SI 2 Cis-inhibiting interactions promote higher 
ratios of high-ligand expressing cells in the Complex 
model. Stationary patterns reached by numerical integration of 
the dynamics for different cis-interactions strengths as measured 
through the cis-binding rates k' h values below the panels (in au -1 
hr -1 units). Ligand levels are represented in grayscale (black for 
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250 and white for 0). Lower cis-binding affinities (k' b =40Q) allow 
high-ligand expressing cells next to each other [25]. Higher cis- 
affinities drive a gradual increase of the ratio of ligand-positive 
cells in the tissue. Herein this phenomenology occurs even in the 
absence of cooperativity (h=l). Parameter values are in the cis- 
inhibition regime: k' s = 0, fi/ = 20 au hr~ 1 , fi n = 5 au hr~ 1 , k s = 25 
auhr 1 , k h = 0A au" 1 hr- 1 , v c = v h = Q.\ hr" 1 , v/ = 0.031 hr" 1 , 
v„ = l hr" 1 , v. 5 = 2 hr" 1 , k- b =k'_ b =0.1 hr" 1 , 6=1000 au" 1 , 
h = 1 . hr refers to hours and au refers to arbitrary concentration 
units. Precursor cells (initial conditions) were set as 
L i (t = 0)=U ! i P,/vi and Ni(t = 0)=U?PJv n where Uf is a 
uniform random number between 0 and 1, and the remaining 
variables were set to 0. 
(TIF) 

Figure S13 Cis-inhibition with a primary Notch signal- 
ing source can create cell-autonomous Instability in the 
Complex model. Representation of relations S9a-S9b in the 
phase space of the signal and the ligand levels. Two stable 
solutions are shown (filled circles) and an unstable solution (empty 
circle). Stability was evaluated through numerical integration of 
the dynamics. This bistability occurs even in the absence of any 
cooperativity (h = 1). Parameter values in the cis-inhibition regime: 
^ 0 = M _ 0 =0.1 hr" 1 , k' s = 0, h = l, fi„ = l au hr" 1 , jg ; = 10 au 
hr- 1 , 6=1000 au" 1 , ^ = 10 au" 1 hr" 1 , k'_ h = 0A hr" 1 , 
Vc = V6 = 0.1 hr" 1 , v„ = l hr" 1 , v s = k s = 0.5 hr" 1 and V/ = 0.031 
hr -1 . hr refers to hours and au refer to arbitrary concentration 
units. 
(TIFF) 

Figure S14 Cell labeling scheme. Arrays of N X M perfect 
hexagonal cells with the subindex labeling schemes used that 
number each cell along the array. In (A) one subindex is used, 
while in (B) we use two subindices. The two main spatial directions 
of the cellular array are depicted. 
(TIFF) 
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